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Abstract 



We study Fredholm determinants related to a family of kernels which describe 
the edge eigenvalue behavior in unitary random matrix models with critical edge 
points. The kernels are natural higher order analogues of the Airy kernel and 
are built out of functions associated with the Painleve I hierarchy. The Fredholm 
determinants related to those kernels are higher order generalizations of the Tracy- 
Widom distribution. We give an explicit expression for the determinants in terms 
Q^l of a distinguished smooth solution to the Painleve II hierarchy. In addition we 

compute large gap asymptotics for the Fredholm determinants. 



1 Introduction 



In unitary random matrix ensembles with a probability measure of the form 

^ ■ ±e-"*'-^WdM, (1.1) 

I on the Hermitian n x n matrices, where V is real analytic on R with sufficient growth 

at infinity, various correlation functions of eigenvalues can be expressed in terms of the 
. kernel 

O ■ Kn{x,y) = {pn{x)pn-l{y) - Pn{y)Pn-l{x)), (1.2) 

X - y Kn 



constructed from the polynomials 



^ . pk{x) = HkX^ H , Kfe > 0, 

orthonormal with respect to the weight e~"^ on M. 

The limiting mean eigenvalue density is known, see e.g. [15], to be given as the 
density of an equilibrium measure ^ly minimizing the logarithmic energy 

/v(/i) = jj log ^^d^(x)ci^(y) + j V{x)dii{x), (1.3) 

among all probability measures ^ on M. The equilibrium density depends on the po- 
tential V but can in general be written in the form [19] 

dx TT ^ " 0, otherwise, 



where Qv is a real analytic function determined by V . Generically Qv has simple zeros 
at the endpoints of supp'i/'V) so that V'v vanishes as a square root [36]. For special (non- 
regular or critical) Vs however, the limiting mean eigenvalue density vanishes faster. 
In general, Qy has a zero of order 4A: + 1, /c = 0, 1, . . . at an endpoint of the support. 
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Remark 1.1 Qv cannot have zeros of order Ak + 3. These endpoint behaviors would 
contradict variational conditions that follow from the minimization property of the 
equilibrium measure [20]. 

Example 1.2 The simplest non-regular case = 1 is realized, e.g., for a critical quartic 
potential 

V(.) = ^^.*-±.' + \.- + l.. (1.4) 

Here the limiting mean eigenvalue density is supported on [—2, 2] and given by 

^^^"^^ = lU;^^^ + ^^'^'^^ " X)^/2X[-2,2](X). (1.5) 

It is easy to verify (1.4) by substituting this i^vix) into the variational conditions. 
In fact, it is much easier to recover V starting from a given limiting mean eigenvalue 
density than to find the limiting mean eigenvalue density corresponding to a given 
potential V directly. 

In order to construct polynomial potentials V for which > 1, it is necessary that 
the degree of V is at least 2k + 2. 

In the generic case where /c = for the rightmost endpoint h = by ol supp'^y, it 
follows from the results of Deift et al [20, 16] that, for any fixed u and v, 

hm ^K„(6 + ^,6 + ^) = if(°)(«,t;), (1.6) 

where 

jX0)(^^ _ Ai (n)Ai^(^) - Ai {v)M'{u) 

' U — V 

is the Airy kernel, and c = cy is a constant depending on V . Let A„ be the largest 
eigenvalue of a random matrix M. It was proved in many cases [39, 16] and is believed 
to hold for = in general that the limiting distribution of A„ is given by 

lim Prob fcn^/2(A„ -h)<s]= det(/ - xf^), (1.8) 

n— >oo \ / 

where K^^ is the Airy-kernel trace-class operator acting on L^(s,(X)). 

The function at the r.h.s. of (1.8) is known as the Tracy- Widom distribution. Tracy 
and Widom [39] discovered a representation 

det(/ - Kf^) = exp j^iy - s)ql{y)dy^ (1.9) 

in terms of the Hastings-McLeod solution of the Painleve II equation 

qxx = xq + 2q^ (1.10) 
characterized by the asymptotic behavior 



go(a^) ~ Ai (x), as x ^ +oo, (l-H) 

1 

8^ 



9o(x) = J ^ ( 1 + ^ + O ) , as X ^ -oo. (1.12) 
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The Tracy- Widom function det(I — K^'^^) does not only describe the largest eigenvalue 
distribution in random matrix ensembles, but appears also in several combinatorial 
models, for example, it provides the distribution of the longest increasing subsequence 
of random permutations [3] . 

In the general case of arbitrary k, a family of limiting kernels K^'''> {u, v) appears in 
place of the Airy kernel. Consider a potential V{x) such that Qv{x) has a zero of order 
4fc + l, = 0, 1, . . . at the rightmost endpoint of suppV'v- Furthermore, to obtain more 
general results, consider for suitable Vj (x) the deformation 

V{x)=V{x)+Y,T,V,{x), 

where Tj are constants. The natural analogues of (1.6) are double-scaling limits where 
n — ^ oo and simultaneously 2} — at an appropriate rate in n characterized by 
parameters tj (in the simplest case we can take all Tj = tj = 0). We refer to [12] 
for more details about those double-scaling limits in the case k = 1. In these limits 
we expect the following expressions for the kernel (corresponding to V{x)) and the 
distribution of the largest eigenvalue, according to conjectures in the physics literature 
[6, 8]: 

for certain constants h and c and any fixed u, v. Here 

—2m{u — V) 

where the functions ^^^''\w) = (^J'^\w\tQ, . . . ,t2it-i) are described below. Moreover, 

lim Prob ( cr?'^^^+^\\n - h) < s\ 
n— >oo \ / 

= i™o^^*(^ ~ ^nX(6+s/[cn2/(4fc+3)],+oo)) = det(/ - K^^), (1.15) 

where KnX(a,b) is the operator with kernel K„ acting on -L^(a, 5), and Kg''^ is the trace- 
class operator with kernel (1.14) acting on L^(s,oo). 

Note that (1.13) for the case k = 1 was proved in [12]. Our goal in this paper is not 
to prove (1.13) and (1.15). Rather, we study the properties of det(/ — i^i^"*). 

The functions ^f^^ = <!>f''\C;to, . . . Mk-i) and •^f'^^ = ^f''\CA), ■ ■ ■ Mk-i) ap- 
pear as solutions of the Lax pair associated with a distinguished solution to the 2A;-th 
member of the Painleve I hierarchy [6, 8, 12]. They are most easily characterized in 
terms of the following Riemann-Hilbert (RH) problem. 

RH problem for $ 

(a) $ = : C \ r ^ C2x2 is analytic, with 



— ITT 



r = U^=irj-U{0}, ri=M+, r3=R-, T2 = e'^iW, T^ = e'^\ 
oriented as in Figure 1. 
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Figure 1: Contour for the $-RH problem. 



(b) $ has boundary values as ( approaches F from the left (w.r.t. the direction 
shown by an arrow), and from the right. They are related by the jump 
conditions 

$+(C) = $-(C)(_°i J), for Cera, (1.16) 

$+(C) = <l>-(C)(j I), forCGTi, (1.17) 

$+(C) = $-(C)(| J), forCGr2UF4. (1.18) 

(c) $ has the following behavior as C — ^ oo: 

$(c) = r ^"^AT (/ + $^c-'/' + o{c-')) 6-^(0-3^ (1.19) 

where $oo is independent of C (its explicit expression will not be important below) , 

'1 0' 
.0 -1, 



'1 

£73 is the Pauli matrix ( ) , is given by 



1 , 



and 

where the fractional powers denote (as usual throughout this paper) the principal 
branches analytic for G C \ (— oo, 0] and positive for ^ > 0. 

(d) $ is bounded near 0. 

The functions $i = and $2 = $2^''^ appearing in (1.14) are the analytic extensions 
of the functions $11 and $21 from the sector in between Fi and F2 to the entire complex 
plane. (One verifies the analyticity by multiplying the jump matrices: see Remark 2.2.) 
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A solution $ to this RH problem satisfies a linear system of the form 



(1.22) 



where ^ is a polynomial in ( of degree 2k + 1, and Bj is a polynomial in (, of degree 
j + 1. The compatibility conditions for the system (1.22) yield the Pp^^ equation, 
the Painlcve I equation of order 4k. The most general form of the RH problem for 

(2k) 

PJ involves a contour consisting of 4/c + 4 rays instead of 4. Our RH problem is the 

(2k) 

one corresponding to one particular solution of Pj , which was introduced by Brezin, 
Marinari, and Parisi in [8]. The precise structure of the Painleve I hierarchy and the 

(2k) 

relation of the functions $i and $2 to the Pj equation is not important for us, but 
we refer the interested reader to [35, 32, 33]. The Frcdholm determinant det(/ — Kg'^^) 
is the object we want to study, and we only need to know the RH characterization of 
$1 and #2 for that purpose. 

Remark 1.3 In the case k = 1, the existence of a RH solution $ was proved in [11, 
Lemma 2.3] for all real values of to,h. The proof of this fact is however valid for 

arbitrary k and arbitrary real values of the parameters tj. It is important to stress that 
the proof does not hold for odd members of the Painleve I hierarchy, in other words for 
the function $(^'^+1). Moreover, the RH problem for is not solvable for every real 
to, because this would imply [21, 11] the existence of a real-valued Painleve I solution 
without poles on the real line, which contradicts [29]. However, as a consequence of 
Remark 1.1, only the kernels generated by the even members of the Painleve I hierarchy 
appear in unitary random matrix ensembles of the form (1.1). 

Remark 1.4 It should be noted that the RH problem for $(2'^) is not uniquely solvable: 



• • • ) t2k-i but not on (. Indeed it is obvious that the left multiplication by a constant 
matrix (independent of does not modify the jump conditions. If the matrix is lower- 
triangular with 1 on both diagonal entries, neither the asymptotic condition (1.19) is 
violated. One verifies as in [28] that, up to the left multiplication by such a matrix, 
the RH problem is uniquely solvable. However, the kernel (1-14) is independent of the 
choice of a solution and thus is well-defined. 

It follows from (1.14) and (1.19) that 



for some constant c > 0. Consequently, the asymptotics of lndet(/ — Ks ) can be 
obtained by a standard series expansion. In particular, 

4fc+3 

lndet(/ - i^f )) = 0{e-'''^), as s +00. (1.24) 

The question of the asymptotic behavior of det(/ — .fci^^) as s —00 is much more 
subtle and will be addressed below. 

1.1 Statement of results 

In the present paper, we will answer the following questions. What are the "large gap" 
asymptotics for the Predholm determinant det(/ — Kg'^'') as s ^ —00? Can we find an 
identity for det(/ — K^^^) which generalizes the Tracy- Widom formula (1.9) to general 
fc? What is the analogue of the Hastings-McLeod solution of Painleve II for general k, 
and what are its properties? 




where lo can depend on 




as u,v ^ +00 



(1.23) 
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1.1.1 Large gap asymptotics 

The first goal of this paper is to find large gap asymptotics for the Fredholm determinant 
det(/ — Kg''^) as s ^ — oo. Those asymptotics are known for the Airy kernel, i.e. in 
the case A; = 0, namely 

lndet(7 - = - 1 In |s| + X + 0{\s\-^/^), (1.25) 

where 

X = -^ln2 + C'(-l), (1.26) 

and C(s) is the Riemann zeta-function. The first two terms at the r.h.s. of (1.25) follow 

easily from (1.9) and (1.12). The expansion of the derivative of (1.25) was obtained by 
Tracy and Widom who also conjectured the value (1-26) for the constant x [39]- A full 
proof of (1.25), (1.26) was given recently in [17] and another proof followed shortly in 
[2]. 

We obtain 

Theorem 1.5 Let K^^^ he the kernel defined in (1-14) for arbitrary k = 0,1,2,..., 

(k) 

depending on parameters to,. . . ,t2k-i, and write Ki for the trace-class operator with 
kernel K^^^ acting on L'^{s, +oo). The asymptotic expansion for the Fredholm determi- 
nant det(/ — Ks^^) as s —>■ —oo is given by the formula 

i >n<iet(/ - K^) = i4(^)W-- + ^ + OIW-*^). (1.27, 



Here ao(s, ■ ■ ■ ) ^2fc-i) and ai(s, • • • > ^2fe-i) are defined as follows: 
aj{s) = 



r(2A; + |) 



r(i + i) \r{2k + 2-j) 

m=2 ^ / 

where the sum vanishes if j = 2k, j = 2k -\- 1, and T(x) is Euler's T -function. 

Remark 1.6 Expanding the right-hand side of (1.27) in negative powers of \s\ and 
integrating, one easily finds asymptotics for det(/ — Kg^^) itself except for the constant 
of integration x^'^^ ■ particular, to the leading order 

lndet(J- kW) = -( TT^ [1 + 0(1)1, (1.29) 

V s ) \2'^k^^{2k)\{2k^\)\) 4fc + 3^ ^ ' 

which was conjectured in [7, 9]. The constant x!"^^ has a representation, as in the Airy 
case A; = 0, in terms of a solution of a Painleve equation and also in terms of a limit 
of multiple integrals (see Section 6). For A; = 0, the relevant multiple integral can be 
reduced to an explicitly computable Selberg integral, which allowed the authors in [17] 
to obtain the simple expression (1.26). 



Example 1.7 For = 0, we have ao(s) = 1 and a\(s) = |, so that (1.27) takes the 
form 



^^lndet{I-KP) = \\sf + ±- + 0{\s\-l). 
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Integrating this expression gives (1.25) without fixing the constant x- 
For k = l, 

5 5 2 

ao(s; to, ti) = - + ti\s\~'^ + 2to\s\~^, ai{s; to, h) = ^ + ^ti\s\~^. 

Theorem 1.5 then gives 

^lndet(I-ii:«) = ^|s|^ + ^ti|s|V^io|s|^ + ^i?|s|' + ioii|s|+tg+gj^ + C(|s|-3), 
which integrates to 

lndet(7 - kW) = -^\s\'- ^t^' - ^to\s\^ - ^tl\sf - hot,\sf - tl\s\ 

+ x(^)-^ln|5| + 0(|s|-2), (1.30) 

where x^^^ is the constant of integration. In Section 6, we discuss 2 representations for 
this constant mentioned in Remark 1.6. 



1.1.2 The Painleve II hierarchy 

On our way to the higher order analogues of the Tracy- Widom formula, let us first de- 
scribe the Painleve II hierarchy and in particular one distinguished solution to it, which 
is a natural higher order analogue of the Hastings-McLeod solution to the Painleve II 
equation. 

(n) 

The n-th member of the Painleve II hierarchy, the PJ^ equation, is the following 

n > 1, (1.31) 



differential equation of order 2n for q = q(x; ri, . . . , r„_i) [38, 14, 30, 35]: 

(d \ ( d \ 

— + 2gj Ln\qx - 9^] + ^ n + C.l\(lx - 9^] =xq-OL, 

where the operator Ln is defined by the Lenard recursion relation 



d_ 
dx 



Cof 



1 

2' 



0, j > 1. (1.32) 



In particular, the first members of the hierarchy (we write qx^ for the j-th derivative 
of q with respect to x) are given by the formulas: 



II 

>{2) 
II 

>(3) 
II 



xq — a, 

Wq^qxx+Gq^) + Ti{q 



xq — a, 
^2 



(1.33) 
(1.34) 



qxx - 2g 

{qi^^ - lOqqi " Wq'qxx + 6(7^) + Ti(fc - 2q' 
(gi^) _ 14(^2^(4) _ 56qqxq^^^ - 7Q{qx?qxx - ^2q{q, 

+ 70q%x + 140g3g2 _ 20g7) + ^2((?i^) " lOgg^ _ iQg^g^^ + Qq^) 

+ Ti{qxx - 2q^) =xq- a. (1.35) 

We are interested in the odd members of the above hierarchy, with the value of the pa- 
rameter a = |. One distinguished solution of the Pj^^ equation attracts our attention. 
It has the following properties. 
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Theorem 1.8 Fix ti, . . . , t„_i G M. For n odd and a = ^, there exists a real solution 

q = q{x; ti, . . . , r„_i) of the P^-^^ equation (1-31) which has no poles for the real values 
of X, and the following asymptotic behavior as x ^ ±00; 

, . 1 , 4n + l ^ 

q{x) = — + 0{x 2^r), as x^ +00, (1.36) 



= |^^^-yy|x|J +0{\x\~^), as x -00. (1.37) 

Remark 1.9 For n = 1, (1.36) and (1.37) together fix the solution to be the Hastings- 

McLeod solution of Painleve II [26] . Note that the function q differs from the function 
u given by (1.10)-(1.12), which is the Hastings-McLeod solution corresponding to the 
value a = 0, while q corresponds to a = ^. We will explain below why the case cx = ^, 
and not a = 0, is relevant for the Predholm determinant if A; > 1. For general odd n, 
although the f[^^ equation is of order 2n, it is our belief that (1.36) and (1.37) together 
still (as for n = 1) determine the real solution q uniquely. We do not have a proof of 
this fact, but we will provide some supporting heuristic arguments later on. 

Remark 1.10 The asymptotics (1.36) at +00 indicate that the terms at the right-hand 
side of (1.31) are dominant compared to those at the left-hand side, and that there is a 
balance between the terms xq and a = ^. There is a whole family of solutions sharing 
the local asymptotic behavior (1.36), see [30]. The asymptotics (1.37) at —00 indicate 
a balance of the g^""''^-term at the left-hand side with the term xq at the right-hand 
side. Again the existence of solutions with this local behavior was shown in [30]. 

Remark 1.11 It is a well-known fact that there is a one-to-one map between the so- 
lutions of a Painleve equation and a so-called monodromy surface [23]. The solution 
considered in Theorem 1.8 is mapped to the point on the monodromy surface corre- 
sponding to the Stokes multipliers 

S2 = S3 = ■ ■ ■ = S2n = 0, 
Si = S2n+l = -i- 

We will characterize this Painleve transcendent in Section 5 in terms of a RH problem. 

1.1.3 Painleve expression for the Fredholm determinants 

Let n = 2k + 1. Define bj = bj{s;to, . . . ,t2k-i) for j = 0,...,2k slightly modified 
compared to aj (1.28) by 

r(2fc + 1) „2fe+i-j 



h.(s) = ^ ' ^' 



2fc-l 



and let x = x{s; to, ... , t2fc-i) and Tj = Tj{s; to, ... , i2fc-i) be polynomials in s defined 
by 

4k+l 

x{s) = -2«+36o(s;to,...,t2fc-i), (1-39) 
Tj{s) = {2j + l)2^^W^bj{s;to,...,t2k-i), j = l,...,2k. (1.40) 

For A; = all ij's, t/s, and the sum in (1.38) vanish. 
We prove 



8 



Theorem 1.12 Let K^^^i he the kernel defined in (1.14) for arbitrary k = 0,1, . . . , and 
write Kg^^ for the trace-class operator with kernel K^^^ acting on L'^{s, +oo). Then the 
following identities hold: 

^ lndet(/ - i^f )) = Q[xisy,nis), T2k{s)], (1.41) 

where x{s),ti{s), . . . ,T2k{s) are given by (1.39)-(1.40). Furthermore, for any fixed 

n, ■ ■ ■,T2k, 

/X 
u{^;Ti,...,T2kfd^ (1.42) 
-oo 

with u{x) = u{x; ti, . . . , T2k) given by 

f + 00 



u{x) 



2-i$^exp|-^ ~(^g(^)_^^dc|.exp|^%(Ode}, (1-43) 



where q{x; ti, . . . , T2k) is the solution of the pj^'^^^^ equation characterized by the Stokes 
multipliers 

S2 = S3 = ■ ■ ■ = S4fe+2 = 0, Si = S4k+3 = 

This solution satisfies the conditions of Theorem 1.8. 

Remark 1.13 By (1.39), the limit x — ±00 corresponds to s ^ ^00. 

Remark 1.14 The function u{x) with fixed t\, . . . ,T2k can also be characterized as 
the solution of the linear second order differential equation 

u"{x) = [q^{x)+q{xf]u{x), (1.44) 

with boundary conditions given by 

4fc+l fx 

u{x) = 2 4fc+3 + 0(1))^ as a; ^+00, (1.45) 

n(x) = 0(1), as a; — *■ — 00. (1-46) 

This allows us to check that our result reproduces the Tracy- Widom formula (1.9) for 
A; = 0. In this case x = — and furthermore, there exists a Backlund transformation 
relating q (solving Painleve II with a = ^) with the usual Hastings- McLeod solution 
go (corresponding to Painleve II (1.10) with a = 0). Indeed we have [4, 14] 

2-^/\x + 2q\x) + 2q,{x)) = q^{-2-^'\f . (1.47) 

Equation (1.44) thus becomes 

u"{x) = [2V3go(_2-V3a.)2 _ |]^(^). (1.48) 

The substitution u{x) = 2^^/^ qQ{—2''^^^ x) = 2~^/^qQ{s) reduces this equation to the 
Painleve II equation (1.10) for qQ{s) w.r.t. the variable s. Moreover, we conclude from 
the boundary conditions (1.45), (1-46) that qo{s) is the Hastings- McLeod solution. Now 
taking one more s-derivative of (1.41) gives 

^\ndei{I-Kf)) = -q^{sf, 
which is equivalent to the Tracy- Widom formula (1.9). 
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Remark 1.15 The function qx{x) + q{x)'^ appearing in (1.44) is a solution of a differ- 
ential equation which is also known as the {2k + l)-th member of the Painleve XXXIV 
hierarchy, and which is equivalent to the {2k + l)-th member of the Painleve II hier- 
archy, see [14] and also [4, 28]. Since the Painleve II hierarchy is more standard and 
better known, we prefer to state our results in terms of it. For k = 0, the function 
Qx + can, as explained above, be expressed in terms of the Hastings-McLeod solution 
Qo for a = 0, relying on the Backlund transformation (1.47) which relates Painleve II 
for a = 1/2 with a = via the Painleve XXXIV equation [14]. For k > 0, a transfor- 
mation for the Painleve II hierarchy relating the values oc = \ and a = is, to the best 
of our knowledge, not present in the literature. 

Outline for the rest of the paper 

The proofs of Theorems 1.5 and 1.12 are based on a RH representation for the logarith- 
mic derivative of the Fredholm determinants det(J — K^^^). We present this identity 

in Section 2. In Section 3, we obtain asymptotics for ^lndet(/ — K^'^) as s ^ — oo 
by applying the steepest descent method of Deift and Zhou to the associated RH prob- 
lem. Section 4 is devoted to the proof of Theorem 1.8, which relies on the study of a 
RH problem for the Painleve II hierarchy. In Section 5, we identify the RH problem 
associated to the Fredholm determinant with the one related to the Painleve II hierar- 
chy. This enables us to give higher order analogues of the Tracy- Widom formula and 
thus to prove Theorem 1.12. In Section 6 we discuss the multiplicative constant in the 
asymptotics of det(/ — ici*^^) as s — — oo. 

2 Fredholm determinant in terms of a RH solution 

In this section, we give a differential identity for the Fredholm determinant det(I — Ks) 
in terms of the solution of a RH problem. This type of identity was obtained in a very 
general framework in [27, 18, 4, 5]. For the convenience of the reader, we recall briefly 
how to do it. 

Recall first that Kg is an integral operator acting on L^(s,-|-oo) with kernel K = 
K^^^ which we can write in the form 



where $i and $2 are related to the Painleve I hierarchy as explained in the introduction. 
Note that 




u — V 




(2.1) 




) 



= {{I-Ks)-\Ks-I + I)){s,s 



) = Rs{s,s), (2.2) 



where I + Rg is the resolvent of the operator Kg given by 



I + Rs = {I- Ks) 



-1 



(2.3) 



The operator Rg has a kernel of the form [18, Lemma 2.8]: 



Rs{u,v) 



{u)H{v) 



F = {I-Ks)-'f, 



H={I-Ks)-^h. 



(2.4) 



u — V 
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Figure 2: Regions I, II, III, and IV. In the case s < 0, at the left, and in the case s > 0, 
at the right. 



The functions F and H can be expressed [4] in terms of the solution 



(2.5) 



of the following RH problem. 

RH problem for Y 

(a) Y is analytic in C \ [s, +oo). 

(b) Y{z) has LP' boundary values related by the condition l+(x) = y_(x)uy(x) for 
X G (s, +oo), with 



(c) Y{C) = I + 0{C-^) asC^oo. 

For n G M, the functions F and H given by (2.4) can be written as [18, Lemma 
2.12] 



where Y{u) is the solution of the above RH problem. 

The functions and $2^^^ now appear in the jump matrix vy- A straightforward 
transformation produces a RH problem with constant jump matrices. Namely, set (see 
Figure 2) 



vy{x) = I — 27rif{x)h (x). 



(2.6) 



F{u) = Y+{u)f{u), 



H{u) = (F-i)^(n)/i(u), 



(2.7) 



^(0^(0, 



for C in region I, II, HI, 




for C, in region V, 



for in region IV, 



if s < 0, (2.8) 



and 



^(0^(0, 



for C in region I, II, HI, 




for C in region IV, 



for C in region V, 



if s > 0. (2.9) 
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Then it is easy to verify that X satisfies the following RH problem [4] . 

RH problem for X 

(a) X is analytic in C \ S*^'*), where 

s(^) = sS^) U s^") U U {s}, 

= (-oo,s), Sj^^ = s + e3F^R~, S|^^ = s + e3^R~, 
oriented as in Figure 2. 

(b) The boundary values of X are related by the jump conditions 

for C G Sl"^ U (2-10) 
forCeS^'^. (2.11) 

(c) X has the following asymptotic behavior as C — ^ oo: 

X{Q = C-i'^^N (/ + XooC"'/' + OiO) e-^(^)-3, (2.12) 

where X^o is independent of ( (its explicit form will not be needed below), N and 
9 are given by 

"=;^(_\ ;)'"'""• 

and 

(d) X has the following behavior near 0: 

X(C)=0(log|C-s|), asC^s. (2.15) 

To obtain (2.12) in the sectors IV and V, note that 

and the real part of 9{Q in those sectors is negative for large \(\. 

Remark 2.1 The existence of a RH solution X is, for real values of s,to, . . . ,t2k-i, 
not an issue, since we constructed X explicitly using Y and Concerning uniqueness 

of a solution X, we note that Remark 1.4 also applies to the RH problem for X, so X 
is only unique up to the left multiplication with a lower-triangular matrix of the form 




x+(c) = x_(c)(J J), 
x+(C) = x_(c)(_°i J), 
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Remark 2.2 Multiplying the three jump matrices for X in counterclockwise direction, 
we obtain 



(2.16) 



(s) 

It follows from this observation that the first column of X in the sector between T,\ 
and (s, +00) has no branching. In other words, the elements in the first column of X 
can be extended from this sector to entire functions in the complex plane. 

Now we will express the logarithmic derivative of det(/ — Kg) in terms of X. By 
(2.7), one verifies that 



^i(c) ; ' 

where Xi and X2 are the analytic extensions of the functions X\i and X21 from the 

sector between T,^^'^ and (s, +co) to the whole complex plane, cf. Remark 2.2. Using 
(2.4) with V = s and u — > s, we then obtain the following identity by (2.2): 

^lndet(/-i^.) = ^ (^-^(C)^'(C)),J,_, (2.17) 

where the limit is taken as C, approaches s from the sector I if s > (sector IV if s < 0) , 
see Figure 2. This differential identity is the starting point of our analysis, and is inde- 
pendent of the chosen RH solution X. In the next section we will obtain asymptotics 
for X i\s s ^ —00, which will lead to asymptotics for the Fredholm determinant. In 
Section 5 we will find a Painleve expression for the right hand side of (2.17). 



3 Large gap asymptotics 

We will now apply the steepest descent method of Deift and Zhou to the RH problem 
for X when s ^ — cx). This approach is very similar to the one used by Kapaev in [32], 
and to techniques used several times in [23]. 

3.1 Rescaled RH problem 

Let us first consider the shifted and rescaled version of X: 

T{C; s,to,..., t2k-i) = X{\s\C + s;s,to,..., i2fe-i). (3.1) 

The shift with s is convenient to 'center' the RH problem at the origin. From the RH 
properties of X, we conclude that T satisfies the following RH problem (see Figure 3). 

RH problem for T 

(a) T is analytic in C \ S, where 

S = El U S2 U S3 U {0}, (3.2) 

—in in 

S2 = (-oo,0), Si = e4fc+3M-, S3 = e4fc+3M-. (3.3) 
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Figure 3: Contour E = Ei U S2 U S3, a circle dUo in which the local parametrix needs 
to be constructed, and the r6gions r2x,ri2?^3- 



(b) T has boundary values which satisfy 

r+(C)=r_(C)(j J), forCGEiUEa, (3.4) 

r+(C)=r_(C)(_°^ J), forCGE2. (3.5) 

(c) As C ^ 00, 

r(c) = (|s|c)"3"3Af (i + i^r^/^XooC"^/^ + o(C"^)) e-^(|^|^+^)^^ (3.6) 

where and 9 arc given by (2.13) and (2.14), and X^x, is the matrix appearing 
in the expansion (2.12) for X at infinity. 

Expanding 9{\s\( + s) as ^ — oo for fixed s < 0, we obtain using (2.14), 

^d^lC + s) = |s|^5(C; s,to,..., hk-i) + kl'^TC'/' + c?(r'/'), (3.7) 

where g is equal to 

2k+l 

5(C;^,«o,...,«2fe-i)= 5^(-i)^+S(«)C'+^ (3-8) 

with the principal branches for the fractional powers, and with the coefficients aj given 

by 

I I r(2fc + |) , r(2fc + |-m) , \ 

J =0,...,2A; + 1. (3.9) 

(the sum vanishes for j = 2k and j = 2fe + 1). The explicit expression for the coefficient 
7 = 7(fe, s, • • • ) *2fe-i) in (3.7) will not be important below. 
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3.2 Normalized RH problem 

The next step is to normalize the RH problem at infinity in a suitable way. We do this 
using the ^-function constructed above. Set 

5(0 = r(C)exp{|s|^5(C)a3}, (3.10) 

so that we have 

RH problem for S 

(a) is analytic in C \ S. 

(b) S has boundary values that satisfy the jump relations 5+ = S^vs on T,, with 
^s(C)=(„ fo, J, forCeSiUSa, (3.11) 



^exp{2|s|^5(C)} 1^ 



^5(0 = y , for C G S2. (3.12) 

(c) S has the following asymptotic behavior as C — oo: 

SiC) = i\s\C)-'^'''N (/ + ^ooC"'/' + OiC')) , (3.13) 
with given by (2.13), and 

Soo{s) = |s|-V2Xoo - 7\s\^as, (3.14) 
where Xoo is the matrix appearing in (2.12). 
Using (3.8)-(3.9), one verifies that for sufficiently large negative s 

Reg{C)<0, asCeSiUSg. (3.15) 

Note that the values of the coefficients tj do not matter for this simple but important 
observation, as long as is sufficiently large. The inequality (3.15) implies that the 
jump matrix vs for S is exponentially close to the identity matrix on S1US3 as s ^ —00. 
For ^ near 0, this uniform convergence breaks down. 

3.3 Outside parametrix 

Ignoring the exponentially small jumps and a neighborhood of 0, we are led to a RH 
problem, which we refer to as the RH problem for the outside parametrix p(°°): 

RH problem for p(°°) 

(a) is analytic in C \ (-00, 0]. 

(b) has boundary values that satisfy the following relation on (— oo,0): 
Pf)(0 = ^'i°°^(C)(_'i J). (3.16) 

(c) As C ^ 00, 

^'^'"^(C) = QslO-'^^'N (/ + 0(C-^/')) . (3.17) 

A solution p(°°) is given by simply removing the error terms in (3.17): 

p(°°\C) = {\s\C)-'^'''N. (3.18) 
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3.4 Local pctrametrix 

Consider a sufficiently small disk Uq (of fixed radius, independent of s) centered at the 
origin. The goal of this section is to construct a function P in Uq which satisfies the 
same jump conditions as S, and matches with the outside parametrix p(°°) for large 
negative s at the boundary dUo of Uq. Thus we have 

RH problem for P 

(a) P is analytic in ?7o \ S. 

(b) P has boundary values that satisfy the following jump relations on C/q n S: 

^-^(^) = ^-(^)(exp{2|.|V,(C)} l)' f-CeC/on(S,UE3), (3.19) 
P+(C) = P-(C)(_°^ J), forCGC/onS2. (3.20) 

(c) As s — — DO, P satisfies the following matching condition with p(°°) at the 
boundary dUo of Uq: 

P(C) = P(°") (C) (I + 0(1)) , for C G dUo. (3.21) 

This local parametrix will approximate S near the origin, while the outside parametrix 
approximates S away from the origin. Wc now construct the local parametrix P ex- 
plicitly using Bessel and Hankel functions. 

3.4.1 Bessel model RH problem 

Inspired by the constructions in e.g. [37, 17], we define the functions Ji, J2, and J3: 

■^'^-^^"2" ^ [-nX'/'H^ {iX^/^) ^AV2JI« (iAV2);' ^"^-^^^ 

where Iq, Kq, and denote the usual modified Bessel functions and Hankel functions 
[1] defined on a plane (rather than universal covering) with the cut along (—00, 0], and 
where wc take principal branches of A^/^, analytic off (— oo,0] and positive for A > 0. 
The functions Ji, J2, and J3 are related to each other by the constant matrices: 



Ji(A) = J2(A) (^J J), for AG C\ (-00,0], 



(3.25) 



J2(A) = J3(A) (^_°^ Jj, forA<0, (3.26) 

J3(A) = Ji(A)(^J fj, for AeC\(-oo,0], (3.27) 
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and have the fohowing large A asymptotics: 



1 + 



1 



-1 
-2i 



-2i 



+ C(A-i) 



0-3 



where 



51 = {\ 

52 = {X 

53 = {X 



-TT + e < arg A < vr — e}, 
-vr + e < arg A < vr}, 
-vr < arg A < vr — e}. 



uniformly as A — > oo in sector Sm, (3.28) 

(3.29) 
(3.30) 
(3.31) 



Using these functions, we can now construct the local parametrix. 

3.4.2 Construction of the local parametrix 

We define the parametrix P in the form 

P(C) = E{C)JM^'+^fiC))ew{\s\'-^9{C)<T3} , for C G ^r, 
where E(Q is an analytic function to be determined later on, 

^1 



(3.32) 



{-77^^ < argC < 77^^} ^ Uo, 



4k + 3 



Ak + 3 



{4^:p37r<argC<7r}nC/o, 

4fe + 2 , 
{-TT < argC < -71-^^} n Uo 



(3.33) 



(3.34) 



4A; + 3 

(see Figure 3), and / is defined as follows: 

(2fe+l 

so that / is a polynomial in ( with 

/(O) = 0, /'(O) = a2(s) > 0. 

Choosing Uq sufficiently small, we see that /(C) is a conformal mapping of Uq which 
maps zero to zero, preserves the real line, and maps flj into 5*^ for j = 1, 2, 3, respec- 
tively. Thus we can use the asymptotic expansion (3.28) and the jump relations (3.25)- 
(3.27) for Jm(|s|^''+V(C))- Namely, using (3.25)-(3.27) and the fact that c/+ = ~g- on 
(—00,0), we deduce from (3.32) that P satisfies the jump conditions (3.19), (3.20) on 

snt/o. 

Finally, define the analytic pre-factor E by 

E{0 = |s|(^+i)'^3^-^V(C)^"^ = I „^ + i-iyajisXn . (3.35) 

From (3.9) wc sec that f{() = f{C,s) tends to some /o(C) independent of s as 
|s| CO. Therefore, for ( G dUo, we can use the asymptotic expansion (3.28) as 
s — > — oo to conclude from (3.35) that 



P(0 = p(°°)(C) 



1 



NI^/(C)V2 



-1 -2i 



-2i 



1 



+ Oi\s\ 



-(4fc+3) 



(3.36) 



It follows that P indeed satisfies the RH conditions of the form proposed at the begin- 
ning of Section 3.4. 
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Figure 4: Reduced system of contours T,r independent of s. 



3.5 Final transformation 

Let us define 



= (_|»|i/2(Ars„Ar-i)2i "^S{()P'-'°>(Cr\ toT(€C\U„, (3.38) 

where 6*00 is the matrix appearing in the asymptotic expansion for S and given by 
(3.14), and N is given by (2.13). Then we have 

RH problem for R 

(a) R is analytic in C \ Hr, with Hr as shown in Figure 4. 

(b) R has boundary values satisfying the jump relations R+ = R-Vr on Sr, with 

vr{0 = PiC)P^^HC)-\ for C G dUo, (3.39) 

VRiO = P^^\OvsiOP^^\0-\ for C G \IJo. (3.40) 

(c) R{0 ^ 7 as C ^ 00. 
Outside Uq, we have for some c > 0, 

4fc+3 

^^ii(C) = ^'^°°^(C)^^5(C)i'^~HC)-' = I + 0(e-'=l^l^), as s ^ -00. 
For C G 5i7o, it follows from (3.36) that in the limit as s — —00 we have 

Z;^(C) =/+|s|-^3^^(^)|5|i-3|<,|-i^+0(|5|-4fc-5/2)^ (3_41) 

where, by (3.36), 

By standard analysis, we see that R also has an asymptotic expansion of the form 
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This expansion is uniform in (, and the l.h.s. is analytic, so that one can differentiate 
these asymptotics with respect to 

Collecting coefficients at the same powers of |s| in the jump relation i?+ = R^vr, 
we obtain: 



Ri,+ iO = Ri,-iO+viiC), 



for C(^dUo, 



(3.44) 



where Ri is analytic in C \ OUq. Since R ^ I as \s\ oo, it follows that Ri tends to 
at infinity. Together with the jump condition (3.44) this normalization enables us to 
compute Ri explicitly: 



^i(C) 



1 



27ri JdUo z-C 



dz 



for C G Uo, 



iRes(^;i(C);0), for C G C \ i7o, 



(3.45) 



where ?^(/) denotes the Taylor part of the Laurent expansion of / near 0. Here we 

used the fact that the pole of vi(Q is simple. 
It is straightforward to check that 



RiiO) 
R[{0) 







ai 
8^ 





8ao " 
*\ 

■it 



(3.46) 
(3.47) 



3.6 Large gap asymptotics 

Using (2.17) and the definition (3.1) of T, we obtain the following identity: 



F{s) := -^lndet(7- kW) = ^^-^ (J-^T^) 



ds 



27ri|s| 
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where C ^ in Sector IV, which by (3.10) leads to a similar identity in terms of the 
normalized RH solution S, 



F{s) = (S-'S'.)^^ 



and by (3.37), 

F{s) 



1 



2m\s\ 



02\ 



+ 



1 



{P-^R-^R!^P) 
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Substituting the parametrix given by (3.32) into this expression, we find 

1 



F{s) = "Q^^^l^l"'^' (J-I J() 



+ 



2ilA^o 2m\s 



+ ——- iJ-^E-^R-^R'cEj) 
27r^ s' ^ ^ ' 
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(3.48) 



To analyze the first term, we use (3.22) and the expansions of Besscl functions near the 
origin (see e.g. [1]) to conclude that 



2iTi 



(3.49) 
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This is the leading order term in the asymptotic expansion for ^ lndet(/ — Kg). Next 
we obtain using (3.35) that 

^ {J-^E-^E'^J)^^\^^^ = Q. (3.50) 



27ri|s 

Computing the last term in (3.48) is slightly more involved, as we need to use the 
asymptotic expansion for R here. We have 



27r s 



Let us now take a closer look at i?(0) ^R'AQ). Using the asymptotic expansion (3.43) 



for R we see that 



i?(0)-ii?^(0) = \s\-''^/^R!^{Q)\s\''^/^\s\-^ +0{\s\-^^-^), 
from which it follows by (3.47) and the Bessel expansions near that 



1 



(J-^E-^R-^R'EJ)^^ = -^^4^ + 0(|s|-^). (3.51) 
^ C ^21 16ao(s)|s| VI I ; \ J 



27ri\s 

Summing up (3.49), (3.50), and (3.51), we obtain 

^.nF,.)4„§„)|,-.J|^.0,W-^). (3.52, 
which proves Theorem 1.5. 

4 A solution to the Painleve II hierarchy 

In this section, we study the RH problem associated with the p|j'^ equation (1.31) for n 
odd. First we pose the RH problem in full generality, with a contour consisting of 4n+2 
rays connecting with infinity. Afterwards we consider one particular pj^^-solution q 
with a = ^, for which the jump contour reduces to four rays. We will then argue why 
q has no poles for real values of (x;ti, . . . ,r„_i), and why it is real. Afterwards we 
will perform a rather standard steepest descent analysis of the RH problem, which is 
similar to the one used in Section 3, in order to find asymptotics for q at +oc. We will 
also obtain asymptotics for q(x) as x ^ —oo. This part requires the construction of a 
y-function and is less straightforward. 

4.1 Lctx system for the Painleve II hierarchy 

The P[^^ equation (1.31) arises as the underlying compatibility condition of a linear 
system. This Lax system was found by Flaschka and Newell in [22] for the Painleve H 
equation, and generalized to the Painleve H hierarchy in [13, 34, 38]. 

The linear system consists of the following equations for 2 x 2 matrices Z = 
Z(C;x,Ti,...,r„_i): 

Z'^ = MZ, Z'^ = BoZ, Z'^.=BjZ, j = l,...,n-l, (4.1) 
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Figure 5: Jumps for ^'(C) if = 1. 



where the logarithmic derivatives with respect to and x have the form 

2n 

' M,r^ — (r.n-o -A- ^ 



M = ^MjC^-xcJ3 + ^c7i, (4.2) 

j=0 



Bo=(-/ I), (4.3) 

and the Bj's are polynomials of degree 2j + 1 in The coefficients Mj are independent 
of (. Except for (4.3) above, the explicit expressions for Mj and Bj in terms of ^, q, 
and Tj will not be needed below. 

The system of Lax equations (4.1) is only compatible if q satisfies the P[j'^ equation. 
Equivalently, in order to preserve the monodromy at infinity for the system Z'^ = MZ 

when varying x, ri, . . . , r„_i, it is necessary that q solves the Pj^^ equation. Using 

canonical solutions to the Lax system for a given Pjj -solution q, one can build a RH 
problem. We do this in the next subsection. 

4.2 The RH problem for the Painleve II hierarchy 

For j = 1, . . . , 4n + 2, let Q,j be the sector 

and denote the boundaries of these sectors 

T, = |c|argC = f^7r| for j = 1, . . . , 4n + 2, (4.5) 



4n + 2 



each of them oriented away from the origin. Set T = U^!!]" Tj. As described in 



[22, 23], for a given solution q to pjj^'' one can choose 4n + 2 canonical solutions Zj 
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_7 = 1, . . . , 4n + 2 to the Lax system (4.1) which have the asymptotic behavior 
Zj{Q = {I + 0(C~^)) e-'^^-''^^"^ , as C ^ oo in the sector 

with 

Prom now on for the rest of the paper we assume that n is odd. Then the function 
defined as 

*(C) = ^i(C) = forCe%, i = l,...,4n + 2 (4.6) 

satisfies the following RH conditions. 

(a) ^'(C) is analytic in C \ T. 

(b) ^'(C) has continuous boundary values satisfying ^'+(C) = *-(C)-(j ^^r C G Tj\{0} 
with 

^^^■+^ = U^-i forj = 0,...,n, (4.7) 

hj = (^Q T)' fori = l,...,n, (4.8) 

hn+i+j = Ij, for j = 1, . . . , 2n + 1, (4.9) 

(c) As C ^ oo, we have 

*(C) = (/ + ^ + 0(r')) e-^®(0'^3^ (4.10) 

with 

e(C; X, n, . . . , r„_0 = 4^^(20^"+^ + E ^^^(20^^'-^^ + xC (4.11) 

Remark 4.1 The asymptotic behavior of the functions is valid in sectors which 
are wider than Oj. The precise shape of the contour T is thus not crucial. It is 
important that Tj is a curve connecting with infinity which hes asymptotically in 
the Stokes sector {( : l^q^fvr < arg^ < 4^^vr}. For simplicity, we have chosen T so 
that it coincides with the anti-Stokes lines on which the leading order term of 6(C) (as 
( — > oo) is purely imaginary. 

Although we suppress this in our notation, ^ depends not only on ^ but also on the 
parameters n, a, x, ri, . . . , T2k (which arc present in the pj"'* equation (1.31)) and on the 
Stokes multipliers (si, . . . , S2n+i)- If, for any set of Stokes multipliers (si, . . . , S2n+i), 
the above RH problem has a solution, it follows from the Lax system for ^ that 

q{x;Ti,.. .,T2k) ■■= 2i*oo,i2(a;,Ti, . . .,T2k) = -2i^'oo,2i(a;, n, . . .,T2k) (4.12) 

solves the Pj"^ equation (1.31). 
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in+2 



Figure 6: Contour T 



However, the RH problem can only be solved if the Stokes multipliers satisfy the 
condition 



For any set of Stokes multipliers lying on the "monodromy surface" given by (4.13), 
the RH problem has a unique solution which is meromorphic in x and ri, . . . ,T2k', it 
determines q{x) by (4.12). The (isolated) values of {x,ti, . . . ,T2k) for which the RH 

(n) 

problem is not solvable correspond to the poles of the Pjj -solution. This map between 
the sets of Stokes multipliers (si, . . . , S2n+i) lying on the monodromy surface and the 
solutions of the pj"^ equation (1.31) is one-to-one. 

4.3 Special solution to the Painleve II hierarchy 

The solution of p|"^ which is of interest to us, corresponds to the value of a = ^ and 
the set of Stokes multipliers 

si = S2n+i = -i, (4.14) 



Because of the trivial jumps on the rays T2, . . . , T2n, our jump contour is now reduced 
to 

T = Ti U T2n+1 U T2n+2 U T4„_|_2. 

Consider the following RH problem. 

RH problem for '^i 

(a) ^'(C) is analytic in C \ T. 

(b) The boundary values of ^{Q on T \ {0} are related by the conditions 



Tr(S'iS'2 . . . <S'2n+i<7i) = — 2isin7ra. 



(4.13) 



S2 = S3 = ■ ■ ■ = S2n = 0- 



(4.15) 





(4.16) 





(4.17) 
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(c) ^ has the following behavior at infinity: 



as C — oo, 



where 

e(C;a;,ri, . . . ,r„ 

(d) Near the origin, 

Ici^ icr 



n-l 



(4.18) 



(2C)2^+^+xC. (4.19) 



*(C) = <^ 



O 



o 



jci5 icr 

ici^ icr 
ici^ icr 



' W 1 



±i 1 



as C ^ 0, C e Si for "+", ( e for 
as C ^ 0, C e ^2, 
as C ^ 0, C e S4. 

(4.20) 



Remark 4.2 The conditions (a,b) are a special case of the conditions (a,b) of the RH 
problem of the previous section, corresponding to our choice of the Stokes multipliers. 
The condition (c) remains the same. For our choice of Stokes multipliers, it was proven 
[22, 23] that ^ can be written near in the form: 



for CeSj, 



with 



As 



-i -1 

1 + ip p 

i -1 

1 — ip p 



A, 



Ai 



-V 

1 p 

i 

1 — ip —i J ^ 



for some p e C, and with E analytic at 0. The condition (4.20) follows from this 
formula. 



Remark 4.3 Using Liouville's theorem, one verifies directly that this RH problem, if 
solvable, has a unique solution. 

From now on, wc write ^ for the solution to this particular case of the RH problem 

(n) 

for Pjj , instead of the more general RH problem with arbitrary Stokes multipliers. 
Starting from this RH problem, we will now prove Theorem 1.8. First, we will prove 
reality and smoothness of q following [10], where this fact was proven in the case n = 1. 
Afterwards we will analyze the RH problem asymptotically using the Deift and Zhou 
steepest descent method in order to obtain asymptotics for q as x ^ ±00. 
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4.4 Smoothness and reality of q 

Let us first prove that the Painleve solution q corresponding to our RH problem is 
real. If ^ solves the RH problem given in Section 4.3, then it is easily verified that 

also (Ji^'(C)o"i and ^)a"i are solutions of the RH problem, with ai = 

Because of the uniqueness we have 

*(C) = <TiW)(^i = <Ti*(-C)<Ti- (4.21) 
Consequently, the matrix ^'oo appearing in (4.18) satisfies 

By (4.12), this imphes the reality of q. 

We will now show that q has no poles for real values of x,ri, . . . ,r„_i. In view 
of (4.12), this is equivalent to the solvability of the RH problem of Section 4.3 for all 
real values of the parameters. A proof of solvability is based on the so-called vanishing 
lemma. This technique was developed in [24, 25] and is explained in detail in [23]. 

Lemma 4.4 (vanishing lemma) Let satisfy the conditions (a), (h), and (d) of 
the RH problem for ^ in Section 4-3, and let the asymptotic condition at infinity be 
replaced by the hom,ogeneous condition 

^o(C)e^®(^^^^ =0(r^), asC^oo. (4.22) 

Then = is the only solution to this RH problem. 

The proof of this vanishing lemma was given in [10, Proposition 2.5] in the case where 
n = 1 and @{(;x) = + x( (for general a > —5)- For arbitrary n and G given by 
(4.11) however, the proof remains exactly the same. It should be noted that the proof 
does not remain valid for any other choice of Stokes multipliers. 

By a standard argument [20, Section 5.3], [23, 31], it follows from Lemma 4.4 that 
the RH problem of Section 4.3 is solvable. 

4.5 Asymptotics for q at +00 for n odd 

We will now perform an asymptotic analysis of the RH problem for the special smooth 

(n) 

solution of Pjj as X ^ +00. This steepest descent analysis shows many similarities 
with the analysis we did in Section 3 in order to obtain large gap asymptotics for the 
Predholm determinant as s — ^ — 00. 

4.5.1 Rescahng and normalization of the RH problem 

Let us define 

T(C) = ^{x^Qexp{i@{x^C)a3}. (4.23) 

First of all, note that the multiplication on the right with this exponential is the most 
obvious way to normalize the RH problem at infinity, see (4.18). Secondly, the rescaling 
with the factor is needed in order to balance the leading order term and the linear 
term of 9 as x — > +00. Indeed we have that 

e{x^C,x,n, . . . ,r„_i) = 

L_(2C)2»+i + C + 0(x-2Tr)') , asx^+00. (4.24) 



2n+l 
X 2n 
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Figure 7: Deformed contour T 



Balancing those terms is crucial for our analysis. The coefficient of C in the above 
expression will turn out to have a direct influence on the leading order asymptotics for 
q{x) as X ^ +00. Let us deform the jump contour T (see Remark 4.1) to a contour 
T as shown in Figure 7, where we let the curves Ti, T2n+i and T2n+2, 'Y^4n+2 coincide 
with the imaginary axis on an interval {—i6,iS) for a small fixed 6 > 0, and we let 
Tj be parallel to Tj away from {~i6,id). Away from the origin, we should choose the 
contour so that it lies in the interior of the region where 



One verifies directly that this is the case for small 5, which is also visible from Figure 
8. The deformation near the origin is not crucial but simplifies the construction of 
a local parametrix later on. The condition (4.25)-(4.26) is essential in order to have 
exponentially decaying jump matrices. We obtain the new T corresponding to T in a 
straightforward way by continuing T analytically through T. Namely, if Ji{C), J2n+iiC)j 
J2n+2{C): Jdn+2{C) j ^-^e the jumps of T on T, i.e. = T_J, and Qj for each j = 
1, 2n + 1, 2n + 2, 4ra + 2 is the region bounded by Tj, Tj, and (0, i5{—iy^^), we set 




(4.25) 



(4.26) 




C € Si, S3, S2, S4 



no = { 




C £ ^2n+l 



(4.27) 



TJ4^n+2, 



We then have 



RH problem for T 



(a) T is analytic in C \ T. 
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Figure 8: Contour plot of Im0(x2n(^) for a; > large and n = 3. The shaded areas 
indicate where lm.@{x'2^Q > 0. The contour T should lie in the shaded region in the 
upper half plane and in the white region in the lower half plane. 



T+iO - r_(c) I _.^2i0{x^o 1 



(b) The boundary values of T on T \ {0} satisfy the jump relations 

for C € 

(4.28) 

for Ce (T2n+2UT4„+2) \H(^,0], 

(4.29) 



r+(c)=r_(c) 



1 _ie-2ie{x^0 
,0 1 



r+(c)=r_(c) 



r+(c)=r_(c) 



,0 



1 



0^ 



force (-^<5,0), 
for C e (0,i(5). 



(4.30) 
(4.31) 



(c) As C — oo, we have 



T{0 = I+^ + OiC^) 



(d) Near the origin, 



ici^ icr 
ici^ icr 



as C ^ 0, ibReC > 0. 



(4.32) 



(4.33) 



It follows directly from (4.24) and (4.25), (4.26) that, as a; — +oo, the jump matrices 
for T tend to the identity matrix exponentially fast on the jump contour T except near 
the origin, say in a disk Uq of radius 6 centered at 0. Moreover, the behaviour of the 
jump matrices for a fixed a; as ^ — > oo ensures that the transformation T — > T preserves 
the condition (c). 
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4.5.2 Local parametrix near the origin 

Wc construct a local parametrix in Uq using Bessel functions. The model RH problem 
we use is different from the one we used in Section 3, although these problems are 
related. 

Bessel model RH problem 

Let us define the functions Mi , M2 as follows: 



Ml (A) = -e-'"'^3/4 - : U / — I21 ' 111 , for A € C \ : 

(4.34) 

3/4 A iHi'\Xe-'-) Hf\\e 



(4.35) 



where H^'^ and H^'^ again denote Hankel functions, with their branch cuts on the 
negative real axis. We then have the relations: 

M2(A) = Ml (A) , for A in the upper half-plane, (4.36) 

Mi(A) = M2(A) , for A in the lower half-plane. (4.37) 

Using the asymptotic properties for the Hankel functions, one verifies that M^ has the 
following asymptotics for A ^ 00: 



Mfc(A) 



as A — >■ 00 in sector Sk, (4.38) 

with 

5i = {A : -TT + e < arg A < TT - e}, (4.39) 

^2 = {A : e < arg A < 27r - e}. (4.40) 

Wc will use Ml and M2 as 'model functions' to construct a local approximation to the 
RH solution T near the origin. 

Construction of the parametrix 

We define the parametrix P as follows: 

P(Q = { Mi(x^/(C))e*®(^^f)'^3, forReC>0,Cei7o, ^^^^^ 
M2(x^/(C))e'®(^^^)'^3, forReC<0, Cet/'o 



where 

/(C) = x-^e(x^O- (4.42) 

Prom (4.24), we have that 

/(O) = 0, /'(O) = 1. (4.43) 
Using the RH conditions for Mi and M2 we obtain (for sufficiently small Uq): 
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Figure 9: Contour 'Er 



RH problem for P 

(a) -P(C) is analytic in Uq \ T. 

(b) P satisfies the following jump relations on C/q H T: 

^+(0 = ^-(C) (_2,,i(.^0 i) ' for C e C/o n [0, ^S], (4.44) 

P+(C) = P-iO (^J -2ie-2i«(^^^) j ^ for C G C/o n [-z,5, 0]. (4.45) 

(c) At the boundary dUo of Uq as x ^ +oo, 

i-ra-Z + ^^f,; -f)+0(x-^), (4.46) 

(d) Near the origin, 
P{0(\ T]=o(\i\l asC^O,±ReC>0. (4.47) 



The local condition (d) follows from the behavior of Hankel functions as A — >^ 0. 

Note that the outside parametrix is simply given by the identity matrix since, away 
from 0, the jump matrices for T tend to / everywhere. 

4.5.3 Final RH problem and asymptotics for q{x) as x ^ +oo 

Let us define 

i?(C) = r(C)P(C)-' forCec/o, (4.48) 

R{0=T{C), forCeC\i7o. (4.49) 

From the behavior of T and P near the origin it follows that R has a removable 
singularity at 0. The function R satisfies a RH problem with jump matrices which are 
uniformly close to the identity matrix. 
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RH problem for R 

(a) R is analytic in C \ S^, with as shown in Figure 9. 

(b) R satisfies the jump relations R+ = R-Vr on Sij, with 

vr{Q = P(C), for C e dUo, 

vr{0 = vt{0, fovCe^R\m- 

(c) R{C) ^ / as C ^ oo. 

Outside Uq, we have already observed that 

2n+l 

^ii(C)=t'T(C) = ^ + C(e-'^^^), asx^+oo. 
For C G 5i7o, we obtain from (4.46) that in the limit as x ^ +00, 



(4.50) 
(4.51) 



MO 



1 



-2i 



(4.52) 



This situation is similar to the one in Section 3.5, and we conclude in exactly the same 
way that 



(4.53) 



X 2n 



with 



i?,(C) = ^Res(.,;0) = i 



i -2i 
i 



for C e C \ Uo- 



Using the fact that T{C) = R{C) for C G C \ C/q and (4.32), we obtain 



1 



~ 8x\2i 
so that by (4.12 



i -2i 
i 



+ 0{x 2Tr 



1 4ri+l 

g(2:;ri,...,r„_i) = 2^ + ^(^~~^ 



as X — +00, 



as X — ^ +00, 



(4.54) 



(4.55) 



which proves part of Theorem 1.8 (equation (1.36). 
For later use in Section 5, consider 



1 



C" 



70-3 



1 1 

-1 1 



—j-irias 



*(iC2)e 4 



1 -2^1^ 

1 



(4.56) 



with the principal branch of the root and the logarithm (the branch cut is (—00, 0)). We 
will need to know the asymptotics of Eq^u = lim^\^o -^'iilC; x,ti, . . . , r^-i) as x ^ +00. 
The ^-limit here is taken so that u = iC,^/'^ tends to zero along (0, +ioo). We have 



e 4 



0,11 



\/2 C\o 



lim C^'^ 



(4.57) 



J 11 



From (4.27), (4.48), and (4.53), we obtain as x ^ +00 



E, 



e 4 



0,11 



V2 C\o 



lim C^/^ 



-1 \){l + 0{^-'^))P-{i^~^C^'^) 



1 
-i 1 



11 



Using the definition of the parametrix (4.41) and the expansions for Hankel functions 
near the origin, we finally find 



-E'0,11 = e 4 -/7rx(l + 0(1)), as X — +00. 



(4.58) 
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Figure 10: Deformed contour T for solutions with 4 non-zero Stokes multipliers 

4.5.4 Asymptotics at +oo for other solutions of the Painleve II hierarchy 

The special solution to P[?^ which we considered, i.e. the one with Stokes multipliers 



is not the only solution which has the asymptotic behavior given by (4.55). For n odd, 
any solution for which the Stokes multipliers satisfy the conditions 



shares this behavior at +oo. These facts can be proved in a very similar way to the 
proof we gave for our special solution q. The main difference in the steepest descent 
analysis is a deformation of the jump contour for T to a contour similar to the one 
shown in Figure 10. If one chooses the curves so that they coincide near as indicated 
in Figure 10, the local paramctrix near can be constructed in exactly the same way as 
in Section 4.5.2. Furthermore, if the contour is chosen so that conditions (4.25), (4.26) 
hold, the jump matrices for R still show uniform decay. 

If some of the 'even' Stokes multipliers S2,S4,. .. are non-zero, it might still be pos- 
sible to construct the local parametrix, but there is no way to deform the jump contour 

in such a way that (4.25) remains valid. In the spirit of the steepest descent method 
for Painleve equations, this means (see [23] for several examples of this procedure) that 
one needs to construct a so-called ^-function, which leads typically to a RH problem 
for which the outside parametrix is not the identity matrix, as it was in our analysis. 
The leading order of the Painleve transcendent will then be determined by this outside 
parametrix. This provides a non-rigorous but strong heuristic argument to believe that 
only the solutions with Stokes multipliers satisfying (4.61)-(4.62) share the asymptotic 
behavior at -|-oo given by (4.55). 

4.6 Asymptotics for q at — oo for n odd 

We will now perform a similar analysis for q{x) as x ^ — oo. The main difference with 
the analysis at +oo is that the normalization of the RH problem needs to be done in a 
less obvious way, using a ^-function. 



■51 = S2n+1 = 

S2 = S3 = ■ ■ ■ = S2n = 0, 



(4.59) 
(4.60) 



Sl + S3 H \- S2n+1 = -2i, 



S2 = S4 = • • • = S2n 



= 



(4.61) 
(4.62) 
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Figure 11: Deformed contour T 

4.6.1 Rescaling of the RH problem and deformation of the jump contour 

Let us first make a rescaling: 

T{0 = ^i\x\^C), (4.63) 

and choose a deformed jump contour T as indicated in Figure 11. We let and Ti 

coincide along [0,(o] and T2n+i and ^2n+2 coincide along [— (^QjO]. Here Co is real and 
positive, and we will determine the precise value of Co later on. Away from [— Co^Co]) 
we take each ray in the contour to make an angle of j^y2 ^ith the real line. It will 
become clear later why this deformation is convenient. We now define a matrix T in 
terms of T, in a similar way as T was defined in (4.27). Namely, in the region bounded 
by Ti, Ti, and (0, Co); as well as in the region bounded by T2n+i, ^2n+i^ and (— Co,0), 
the function T is defined by extending T analytically from 82- In the region bounded 
by T4„+2, T^4„+2, and Jo, Co), as weU as in the region bounded by T2„+2, T2n+2, and 
(_^Co;0), the function T is defined by extending T analytically from 84,. In the regions 
^i, 53, ^2 and ^4, we set f = T. 

4.6.2 Construction of the g'-function and normalization of the RH problem 

In the previous section where we computed asymptotics at +00, we used the obvious 

way to normalize the RH problem, namely, multiplying on the right with e*®(l^l^^)°'3. 
Here we could do this as well, but it would not lead to a RH problem with decaying 
jump matrices. This is a consequence of the fact that the topology of the set {C : 
Im9(|x|2;^C) > 0} for negative x is different from that in Figure 8. We will deal with 
this problem by replacing by a function g behaving like G(|x|2;^C) at infinity and 
such that the set {C : Im g{Q > 0} has a convenient topology. 

We define g in the form 

n 

5(C) = E-.(C'-Co^)'^> (4-64) 

so that g{() is analytic in C \ [— Co, Co]- 

We fix the constants Cj and the branch point Co > by the requirement that 

\x\'-^g{0 = e{\x\^0+O{^), asC^oo. (4.65) 
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Expanding this expression as C ^ oo and equating the coefficients at ^^n+i^ ^2n-i^ • • • , C'^ 
gives (after some calculations which make use of the binomial formula and induction) 
that, as X — ^ — oo, 

22n 

r(jj -L. I) o2n-l 

cn-m = ^, ^ ^ , Co"' + 0(1^1""), m = 1, . . . , n - 1. (4.67) 

r(n — m+§) ml 

The coefficient at C yields 



(2n = '-v- + 0(1x1-^) = + 0(1x1-^). (4.68) 
^° 22nr(n + i) ^' ' ^ (2n)! I ^ ^ ' 



Here the Odx] «) can also be computed explicitly by taking into account the terms 
with Tj in 6(|x| 2;^(^). Note that the rescaling with a factor |x|2^ in (4.63) was necessary 
to have branch points ±(o with a non-zero finite limit as x — — oo. 

We now establish some useful properties of the function g{Q. First, let ( G (0, Co)- 
Then 

IC'-Co'l = Co'^, zG(o,i). 

Note that by (4.66), (4.67), and (4.68), 

Cn-mCo - Ti^t , 3x , + C(fI ")> m = 0, ...,n-l, x^-oo. 

2r(n — m + §)m! 

Thus 



^ ±i7r(2j+l)/2 j+l/2 

..(0 = ^C.E r(,- + 3/2)(„-,)! +"«^|-^'' <^-«=" 

Therefore, 

5+(C)-5-(C) = 2g+{0 =iV^nKoct>{z)+0{\x\-Ti), ( G (0,Co), ^ e (0,1), (4.70) 
where 

y tlM . (4.71) 

^r(i + 3/2)(n-i)! ^ ^ 



We need to determine the sign of 0(z) for z € (0, 1). We do it by identifying 0(z) with 
a Jacobi polynomial. Jacobi polynomials p^^'^\x) depend on 2 parameters a and 
and have several representations in terms of the hypergeometric function F. We will 
make use of the following 2 of them: 

^-"'^"^ = in -mill a/ (- + « + A -n + 1; 1 + a; , (4.72) 

^ r(n + /3)(x-l)-^ f 1 - n - a -n + 1- 1 + /3- (4 73) 

2-i(n-l)!r(l + /?) V ''+''' + ^'x-ly' ■ ^^"^^ 

Prom the first of these equations we obtain for q: = 3/2, /3 = —n — 1/2: 

^(3A-n-i/2)^^ _ 2.) = -I>±M^(,). (4.74) 
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On the other hand, the second equation for this particular Jacobi polynomial yields 
f3/2,-n-i/2),, , . _ 2r(n + 3/2) z"— 1(1 - zf 



(3/2,-n-l/2) . _ . _ Xn^oi^^ - . . 

^"-1 ^ ^ m!(n-l-m)!(n + l/2-m)' ^^'^^^ 



m=0 

which is obviously positive for z G (0, 1). Therefore, it follows by (4.74) that 

</>(z) < 0, z G (0, 1), (4.76) 

and hence, from (4.70), 

Im (5+ - 5_) = 2Im5+ < 0, C e (0, Co), (4.77) 

for sufficiently large |a:;| depending on r^'s. 

If C G (— Co,0) the only difference is that and "— " sides are interchanged and 
therefore 

Im (5+ - 5_) = 2Im5+ > 0, C e (-Co, 0). (4.78) 
Moreover, considering the arguments of C — Co and C + Co, we easily conclude that 

Imff(C)>0, CgTiUT2„+i (4.79) 

Im5(C) < 0, C G T2„+2 U T4n+2 (4.80) 

5+(C) + 5-(C) = 0, Ce(-Co,Co)- (4.81) 

Note that the inequalities here and in (4.78) hold, as in (4.77), for sufficiently large 
Let us now define 

S(C) = f (C) exp \i\xf^ g{Qa^^ . (4.82) 
Prom the ^'-RH problem we then have 

RH problem for S 

(a) S is analytic in C \ T. 

(b) S satisfies the following jump relations on T \ {0}: 

«+(« = ^-K'(-iexp{2iN^s(C)} l)' f-<<^(*'uM\l-^»-f°l' 

(4.83) 

S^O = S.iC) (^J -^^^P {-2^N^.(C)} j ^ ^ , (Y^^^^ ^ Y^^^^) \ 



(4.84) 



exp|2i|a;| 2n ^_|_(C)| —i 



S+(0 = S_«)("*-l"'-^, »^>^'J J, for(e(-Co,0), (4.85) 

S.(C) = S-(C)(° for<e(0,<„), (4.86) 
(c) At infinity, 

S{0=I + O{C-')), C^oo. (4.87) 
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(d) Near the origin, 



o 
o 



|C|V2 |^|-l/2\ 
ICI-V2 |^|l/2l' 



as C ^ 0, ImC > 0, 
as C ^ 0, ImC < 0. 



(4.88) 



It follows from (4.77)-(4.80) that the jump matrices for S tend exponentially fast to 
the identity matrix away from the interval [— CojCo]- On (— CoiCo)) the diagonal entries 
decay exponentially fast. The condition (c) is also easy to verify. 

4.6.3 Parametrix outside ±^0 

Ignoring exponentially small jumps and small neighborhoods of the branch points ±^0) 
we are led to a problem for a parametrix 

RH problem for p(°°) 

(a) P(°°) is analytic in C \ [-Co, Co]- 

(b) satisfies the following jump relations on (—Co, 0) U (0, Co), with the intervals 
(— Co,0) and (0, Co) both oriented away from the origin, 



-i 



(c) At infinity, 

pM(C) = / + o(C-'), C^oo. 



(4.89) 



(4.90) 



(d) Near the origin. 



p(oo)(^) 



o 

O 



|^|l/2 |^|-l/2> 
^|C|-V2 |^|l/2^ 



as C ^ 0, Im C > 0, 
as C ^ 0, ImC < 0. 



(4.91) 



It is easily seen that this problem has a unique solution which can be found as 
follows. First, let P(C) satisfy the same conditions (a), (b), and (c) as P(C), and write 
it in the form 



1 fl -1 



2 Vl 1 



1 1 
-1 1 



(4.92) 



-i 



Then L(C) has a jump on (—Co, Co) with the diagonal jump matrix ^ J , and 

L(oo) = I. Therefore, considering the 2 scalar RH problems, we immediately obtain a 
solution 



m = a{Cr\ a(C) 



1/2 



c 



(4.93) 
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We assume that the roots are positive on the real axes to the right of the respective 
branch points, and the cuts go along the real axis to the left. Note that the correspond- 
ing P(C) satisfies the conditions (a), (b), and (c), but not (d). Let us multiply it on 
the left by the following matrix 

and choose a,P,j,S so that AP{() satisfies (d), and therefore, since the conditions 
(a,b,c) are not violated, gives the solution p(°°) we are looking for. A straightforward 
calculation yields a = (3 = —1(0 /2, j = 6 = i(o/2. Thus, we obtain 



As — 00, we clearly have 




(4.95) 



^'^°°^(C) = /+^(/ (4.96) 

The off-diagonal entries in the coefficient at in this expansion will determine the 
leading order asymptotics for q{x). 

4.6.4 Local parametrices near —(0 and ^0 

In the neighborhoods U of Co and U of — Coi the local parametrices are constructed in 
terms of Airy function. The analysis is standard, and we only present the result. In U, 
we have 

' , , /, , 2n+l \ r , ,2n + l , , ~1 

P(C) = ^(0^ (kl^Aj exp {i\x\^giC)a3\ , 

^(C) = P^-HO^: (I 7) (|x|^a)^^^\ A(C)^/^ = le^'-/'giC), (4.97) 

where the minus sign is taken in the upper half-plane, and plus, in the lower. The 
function A has the following form in the sector between T2n+i and the real-axis-part 
of the contour: 

and is constructed in the other sectors using the jump conditions. 

We have the crucial properties: P{()S{()~^ is analytic in U, and P{()P^°°\()''^ = 
I + 0{\x\ 2;^) uniformly on dU as 00. The parametrix P{Q for ^ G J7 is 

obtained using the transformation aiA{ — z)ai. 

Note that these parametrices do not contribute to the leading-order asymptotics for 
q{x). 

4.6.5 Final RH problem and asymptotics for q at —00 

As before we define R by 

i?(C) = SiOPiC)-^ ioiCeUUU, (4.99) 

i?(C) = 5(C)P(°°)(C)-\ forCeC\(C/UC/), (4.100) 

so that R satisfies the following RH problem. 
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Figure 12: Contour for R. 



RH problem for R 

(a) R is analytic in C \ E^, with as shown in Figure 12. 

(b) R satisfies the jump relations R+ = R-Vr on E^j, with 

vr{Q = P{QP^^\0-\ for C e 5C/ U at^, (4.101) 



vr{Q = P^'^\Ovs{OP^'^\0-\ for C G Ei^ \ (t/ U U [-Co + e, Co - e]), (4.102) 
.i.(C) = / + Pi°°^(C)f° --p{2^N'^^+(0}Vi°°^(C)-\ forCG(-Co + .,0), 







(4.103) 



iexp|-2i|x|^£?+(C)| 

(4.104) 



t;;^(C) = / + Pl~)(C) r o.u,2^w^^\ n^-°°^(Cr\ forCe(o,Co-^). 



(c) R{Q ^ / as C ^ oo, 

where £ is the radius of the neighborhoods U and [/. Note that, beeause of (4.88) and 
(4.91), R is bounded near the origin. From the last sections we recall that vr = I + 

. . I _ 2n-|-l > „ „ / N / I 2n+ 1 , 

0(\x\ 2n ) uniformly on Er. It follows in a standard way that i?(C) = -f + 0(|a;| 2n ) 
uniformly in C as x — — oo. Therefore, we obtain from (4.100), (4.96), (4.82), and (4.10) 
that 

*oo = ^W^("/ "/)+0(W-i), (4.105) 
so that, by (4.12), 

g(x) = Cokl^ + 0(|x|-^), asx^-oo, (4.106) 

which, in combination with (4.68), proves the final part of Theorem 1.8. 

For later use, we again need asymptotics for E defined by (4.56). In Section 5 we 
show that E{C,) is analytic near the origin, and we can write 

£;(C) = £;o(/ + £^iC + o(C)). 

A similar analysis to the one following (4.56) for x — > +oo gives 

-^0,11 = o(l), as a; ^ -oo. (4.107) 

and 

^1,21 = o(l), as x ^ -oo. (4.108) 
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4.6.6 Asymptotics at — cxd for other solutions of the Painleve II hierarchy 

Similarly to the situation at +00, the asymptotics at —00 only do not uniquely fix our 

(n) 

special sohition of Pjj . In a similar way as was shown for the Painleve II equation 
itself (see [23]), one has the freedom to choose some of the even Stokes multipliers 
S2, S4, • • • , S2n different from zero without destroying the above RH analysis. 

Under the assumption that the even Stokes multipliers are zero, the method works 
only if the jump contour connects the branch points ±^0 with infinity within every 
Stokes sector corresponding to a non-zero Stokes multiplier. A simple analysis of g(Q 
shows that this only leads to decaying jump matrices for the Stokes multipliers si and 
S2n-i-i) so the other odd multipliers should be zero. Recalling Section 4.5.4, we can 
conjecture that the restriction 

S2 = S3 = ■ ■ ■ = S2n = 

is needed in order to connect the asymptotics of q{x) at +00 and —00. In order 
to preserve the symmetry relation (4.21) and thus to have a real solution, the only 
possibility then is that si = S2n+i = —i- 

Let us stress once again that this is not a proof that the Painleve solution wc discuss 
is uniquely determined by the reality and by its asymptotics. In principle a different 
asymptotic analysis of the RH problem could still lead to the same asymptotics for q. 
However in view of the direct dependence of the leading order asymptotics for q on 
the outside paramctrix, it would be rather surprising if other solutions with the same 
asymptotics existed. 

5 Higher order Painleve II formula for the Fredholm de- 
terminant 

In this section, we will relate the X-RH problem for the Fredholm determinant to the 
^'-RH problem for the Painleve II hierarchy, namely, to the one corresponding to the 
special solution of pj"'' we described in the previous section. By (2.17), this will enable 

us to express the logarithmic derivative of det(/ — Kg'^^) in terms of the Painleve II 
hierarchy and to prove Theorem 1.12. 

5.1 RH problem for the Painleve XXXIV hierarchy and the Fredholm 

determinant 

Let us define the following function U = U^"'^ {^;x,Ti, . . . , r„_i) in terms of the solution 
^'("^ of the RH problem stated in Section 4.3: 

C/(-)(0 = ^r^"^ J) v&(-)«^;x,ri,...,r„_i)e-^-3. (5.1) 

Here we take as usual the branches of the fractional powers which are analytic in 

C \ (—00, 0] and positive on (0, +00). 

To obtain a RH problem for U from that for ^, note first that the ^'-RH problem is 
defined in the plane of the variable u = iC^^^- Because of the symmetry ai^{—u)ai = 
'^{u), it is enough to consider only the half-plane Im-u > 0. The function ^ = {—iuY 
maps it onto the whole C-plane. Now it is easy to verify that U satisfies the following 
conditions. 
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RH problem for U 

(a) U {Q is analytic in C \ S, where the contour S is the one shown in Figure 3, 

S = Si U S2 U S3 U {0}. (5.2) 

(b) The boundary values of U satisfy the relations: 

[/+(0 = C/-(C)(| J), forCeSiUSa, (5.3) 

i7+(C) = C/-(C)(_°i J), forCGS2. (5.4) 

(c) As C — 00, 

U{Q = C^^'N + C»(r^)) e-^®(*«^)'^3^ (5.5) 

with N given by (2.13) and 9 by (4.11). 

(d) As C ^ 0, 

U{0 = O{\C\-'/'). (5.6) 

Although we do not use this fact, it is worth noting that the RH problem for U is 
the RH problem associated to one solution of the 2n-th order equation in the Painleve 
XXXIV hierarchy 

Our aim is to identify U = U^"'^ with the RH problem for X^^^ from Section 2 
if n = 2A; + 1. Let us first identify the parameters {s,to, - ■ ■ ,i2fe-i) for X with the 
parameters {x,ti, . . . ,T2k) for U. We do this by requiring that 

1 /n 1 4k+l 

e{C + s;to,...,t2k-i)=ie{i{cO'^^) + 0{-), asC^oo, = 4-^^, (5.7) 

where 9 is defined in (2.14), and 9 in (4.19). The value of c follows from comparison 
of the leading order terms in ( for 6 and 9. Expanding 9{(^ + s) for — ^ 00 similarly 
as in (3.7)-(3.9) for fixed s, we find 

2k+l 

9iC + s)=Y.bjC'^+0{C'% (5.8) 



with 



r(2A; + 

r(j + |)r(2fe + 2-j)" 



2k-l 



Comparing the same powers of ^ in (5-7) we obtain that x = x{s;to, . . . ,t2fe-i) smd 
Tj = Tj{s; to, ... , t2k-i) are given as follows: 

X{s) = -C-'/Ho{s; to,..., t2k-l), (5.10) 

Tjis) = (2j + l)A-^c-^-h,{s; to, . . .,t2k-i), j = l,...,2k. (5.11) 
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Let us now set 



X(C; to,..., t2k-i) = c4°'3?7(c(C - s);x,Ti,..., T2k), 



(5.12) 



with X and tj given by (5.10) and (5.11). It is easy to verify tliat X satisfies the jump 
conditions (2.10)-(2.11) of the RH problem for X, and therefore XX~^ is a holomorphic 
function in C \ {0} for any RH solution X. Using (5.7) and the asymptotic behavior of 
U, we see that also the asymptotic condition (2.12) holds for X and XX~^ is bounded 
at infinity. By (2.15) and (5.6), XX~^ can only have a removable singularity at 0. 
Therefore, by Liouville's theorem, XX~^ is a constant matrix, and so X solves the RH 
problem for X. Since (2.17) holds for any solution of the X-RH problem, we obtain 



4fc+l 

F{s) = - lndet(7 - i^i^)) = {U-\z)U,{z)) 



ds 



2m 



2iU\o • 



(5.13) 



5.2 Proof of Theorem 1.12 

In order to prove Theorem 1.12, we will consider the x-derivative of U. As already 
mentioned, the RH problem for U (see Section 5.1) is not uniquely solvable but has a 
family of solutions of the form 



where U is defined by (5.1), and rj is independent of ^ but may depend on x, ri. 
Although the identity 



4fc+l 

F(.) = -lndet(/-i^<')) = — 
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(5.14) 



(5.15) 



which we obtained in (5.13), is independent of the choice of r), the linear system satisfied 
by U does depend on rj. A convenient choice for us is 77 = —q. Then we obtain from 
(5.1) and (4.1), (4.3) that U satisfies the following differential equation with respect to 
X (which is one of the equations in the standard Lax pair for Painleve XXXIV): 



UAC;x) = 







-1 



U{C-x). 



< - {qx + q^) 

Let us define a function E by the equation 

'1 h^^c 







1 



c{0. 



(5.16) 



(5.17) 



where 




in region I, 
in region II, 

in region III. 



(5.18) 



The sectors I, II, and III are those indicated in Figure 2 if s = 0. Considering the RH 
problem for U and choosing the branch of the logarithm with the cut along (— oo,0). 
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it is easy to verify that (5.17) defines an analytic E{() in a neighborhood of zero with 
det^ = 1. 

Writing a Taylor series for E near 0, 

E{C;x) = Eo{l + EiC + 0{(:^)), asC^O, (5.19) 

we obtain from (5.15) that 

^lndet(/-iCf)) = g(x;ri,...,r2fe), (5.20) 

with 

4k+l 

Q{x;Ti,...,T2k) = Ei^2l{x;Tl,...,T2k)- (5.21) 

Using (5.16) we also find 

3 0' 
-1 



El,. = E^ ' j ° ^'j Eo, (5.23) 



and thus 



£^1,21,. = -£^0,11- (5-24) 
In other words, 



o 6 4 4fe+l 

Q'{x) = u{xf, u{x) = -=2-«+3£;o,ii(x). (5.25) 

V 27r 



We see from (4.108) that £"1,21 (-oo) = 0. Therefore equations (5.21) and (5.25) 
imply the formula (1.42). 

Using (5.22) we obtain the following second-order differential equation for u{x): 



u 



"{x) = [q.{x)+q\x)]u{x). (5.26) 



Substituting the asymptotics (4.58) and (4.107) into (5.25), we obtain the boundary 
conditions (1.45)-(1.46) for u. They determine u{x) uniquely for a given q{x). Now it 
is easy to check directly using Theorem 1.8 that u{x) given by (1.43) is the solution 
of (5.26) satisfying these boundary conditions. This completes the proof of Theorem 
1.12. 



6 The constant problem 

In this last section we discuss the multiplicative constant in the asymptotics of det(/ — 
as s —00. In contrast to the sine, Airv and Bessel kernel determinants, we 
do not provide a compact expression for this constant in terms of single integrals of 
elementary functions. For this reason and for simplicity of the argument, we assume 
that k = 1 and to = ti = for the rest of this section (a generalization being straight- 
forward). In this case, the expression (1.30) has the form 

lndet(7 - kP) = -^\s\' + x^'^ " ^ In \s\ + 0{\s\-^), s ^ -00, (6.1) 
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where x^'^'^ is the logarithm of the constant in question. Comparing this expression with 
Theorem 1.12 for A; = 1, we obtain a representation of the constant in terms of our 
special solution q{x) of Pjj^ 



X^^^ = li™ { ^ ^ l^r + — In \s 



s— >— oo 



28-7' 

+00 



/+00 p -I \ 

Q [-25/7-35^3; 2^7-215^2^ 2-^/'bt\ dtj , (6.2) 



which is basically (using (1.42) and (1.43)) a closed expression for x^^^ in terms of the 
Painleve II solution q with A; = 1. 

This is a generalization of a similar formula, expressing x^^^ in terms of Pjj\ in the 
case of the Airy-kernel determinant. 

Let us now try to generalize the argument in [17] which led to the simple expression 
(1.26) for x^^^ ■ Consider the kernel (1.2) with V(x) given by 

V{x) = 4 - ^x^ + 3x^ - 2x 

for the polynomials Pk{x) orthogonal with the weight cr^^i^) on the half-axis (0, -|-oo) 
(we choose a half-axis instead of M, and hence this V instead of (1.4), for convenience 
only). The corresponding mean eigenvalue density is supported on [0, 2] with the k = 1 
type singularity at x = 2: 

V'y(x) = |-xV2(2-x)^/2;^[0,2]. 

Remark 6.1 Note that, on the half-axis, one can find a potential V of degree 3 with 
a singular endpoint. However, this does not simplify the argument below. 

We now expect equation (1.15) to hold at the endpoint 6 = 2. To prove the conver- 
gence of the determinants (1.15), one would need to obtain global estimates for the 
polynomials Pk{x) (cf. [17, 16]), which should be possible in view of [12]. Note that 
(1.15) for our V{x) = V{x) is equivalent to the following: 

■im °"P+;/'°;'^'» =det(J-gW). (6.3) 

n^oo Dn{+00) \ s h \ J 

where 

Dn{a) = 1- r... r H {x,- xjf n e-^^'^'-^^^dx,. (6.4) 

In view of (6.1), this formula would provide a representation for x^^^ in terms of a 
limit of multiple integrals. It can be however simplified following [17]. Namely, there 
exists a differential identity for ^ In Dn (a) in terms of the polynomials orthogonal with 

the weight e"^^^^) on the interval (0, a). Analysing the related RH problem, one may be 
able to obtain a uniform asymptotics for the polynomials, and hence, for ^lnD„(a) 
in the range a G (0, 2 -|- ^/(cn^/^)), s < sq < for a sufficiently large [sq], n > \so\. On 
the other hand, it is easy to obtain a series expansion for Dn{a) for n fixed and a ^ 0. 
Combining these results, one could integrate ^ lnZ)„(a) from a ^ to 2 -|- sl{cr?l'^\ 
and therefore one may be able to obtain, as in [17], an expansion for L>„(2-Fs/(cn2/'^)) 
as n — > oo (and \s\ sufficiently large) in terms of elementary functions. Thus, the only 
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non-elementary object to enter the expression for det(/ — KI ), s < sq, and hence the 
expression for x^^\ would be Dn{+oo). This would appear to be a final formula for x^^^- 
Indeed, note that in [17], for x^^\ the corresponding D„(-|-oo) was a Selberg integral 
and therefore the corresponding formulas simplified to (1.26). In the present case, 
however, there is no known expression for D„(-|-oo) in terms of elementary functions 
(or a fixed number of integrals thereof). An attempt to derive the asymptotics of 
Dn{+ oo) as n — oo by continuing the integration of the differential identity beyond 
the endpoint 2 of the measure is likely to produce an expansion involving integrals of 
Painleve functions. 
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